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Abstract 



We present preliminary results on the branching fractions B(B° — ► D* + D*~) and B(Df — > </>7r + ), 
based on a data sample of approximately 124 x 10 6 BB events collected by the BABAR detector 
at the PEP-II e + e~ B-factory. B(B° — ► D* + D*~) is measured selecting neutral B meson decays 
to the final state D*~ D* s + with partial reconstruction of the D* + , in which only the D*~ and 
the soft photons from the decay D* + — * Dj*y are reconstructed. The branching fraction product 
B(B° — ► D* s + D*~) ■ B(D+ — > (/>7T + ) is measured via a complete reconstruction of the whole decay 
chain. Comparing these two measurements provides a model-independent determination of the 
B{Df — > 4>tt + ) branching fraction. We obtain the following preliminary results: 

?o 



B(B° -> D*+D*~) = (1.85 ± 0.09 ± 0.16)% 







B(D+ -> 07T+) = (4.71 ± 0.47 ± 0.35)% 
where the first error in each measurement is statistical, the second systematic. 
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1 INTRODUCTION 



A precise measurement of the branching fraction for the Df — > c/>7r + mode is important because 
nearly all Df branching fractions are determined by normalizing the measurements to B{Df — > 
07r + ) 1 . The present uncertainty of about 25% on B{Df — * 2^ thus affects many of the 

results regarding mesons. 

In the factorization model for two-body decay rates, it is assumed that the transition amplitude 
of the process is the product of two currents that can be evaluated separately. This model has been 
successful (3 in describing the measured branching fractions and polarizations for B meson decays 
such as B° — ► D*~tt + 0], B° — > D*~ p + and B° — ► D*~af t 5j, in which the momentum transfer in 
the process is low (q 2 ~ m^m 2 ). Measurements of decay rates for modes such as B° — > D* + D*~ 
allow tests of the predictions for high q 2 . 



2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring 
and correspond to an integrated luminosity of 112.3 fb -1 and to (124 ± 1) x 10 6 BB pairs. A 
detailed description of the detector can be found in Ref. 0. 

In addition to this data sample, several simulated event samples were used in order to study 
efficiencies and backgrounds. For background studies, we used a Monte Carlo simulation of B®B® 
and B + B~ events (each equivalent to an integrated luminosity of 440fb -1 ), e + e~ — > cc (230 fb -1 ) 
and e + e~ — > uu, dd, ss (180 fb -1 ). 



3 ANALYSIS METHOD 

The aim of this analysis is to measure the branching fractions B(B° — » D* s + D*~) and B(Df — > 
4>tt + ). To accomplish this, the B° — > D* + D*~ (Dfj)(D°ir~) decay is reconstructed using two 
different methods. 

The first method uses a partial reconstruction technique, in which only the D*~ is fully recon- 
structed, and then combined with the soft photon from the D* + — > Dfj decay, without requir- 
ing explicit Df reconstruction. Denoting the measured event yield by Md s i we can express the 
B° — ► D* + D*~ branching fraction as: 



Bi = B{B° -► D* + D*-) = K Ds n0 ■ (1) 

y,(~,-Br< 



Here K = [2N ^)B(D*+ D+-y)B(D*~ £>°7r-)] _1 > iV _o is the number of neutral 5 meson 

pairs, B/^ are the branching fractions for the D° decay mode i, £{ are the efficiencies for recon- 
structing and selecting the partially reconstructed B° into a final state containing a photon, a soft 
pion and a D reconstructed into mode i. We assume B(T(45) — ► B B°) = 0.5. 

B° — > D* s + D*~ decays can also be fully reconstructed. In this paper we focus on the Df — ► (/>7r + 
mode, where the <p meson is reconstructed in the K + K~ channel, determining from the measured 
yield ND a ^<]>-K the product of branching fractions B(B° — > D* + D*~) ■ B(Df — > 07r + ): 



B 2 = B(£° -» d;+D*-)B(£)+ -> <M+) = /C -A/g^frr =5-, (2) 

B(</>^ A+A-)£^-Bf° N 
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where is the efficiency for detecting the fully reconstructed B°, including reconstruction of 
<j> — > K + K~ . The branching fraction B(Df — > 07r + ) is measured from the ratio of the two yields: 

{ s } ~ Bi " M Ds B(4> - *+jFj • W) ' 

where the factor /C exactly drops off, and although the efficiencies E{ and £[ are in general different, 
many systematic uncertainties cancel in the ratio, as will be discussed in Sec. El 



4 PARTIAL RECONSTRUCTION ANALYSIS 

4.1 Signal Extraction 

We reconstruct the B° — > D* + D*~ — > (Df^)(D°ir^) decay by combining photons in the event with 
fully reconstructed Z)* - mesons, without requiring reconstruction of the from the D* + decay. 
In order to extract the signal, we compute the missing mass m m ; ss recoiling against the D*~ — 7 
system (all quantities defined in the center-of-mass (CM) frame): 



m miss = (£beam ~ E D * - E^) 2 - (p B - p D * ~ P 7 ) 2 . (4) 

The 

m miss distribution of signal events peaks at the nominal mass [J with a spread of about 
15MeV/c 2 . The kinematics of the event are not fully constrained with this procedure and one of 
the decay parameters must be chosen in an arbitrary way. In particular, taking the beam energy 
in the CM to be the B energy, the angle between the B momentum vector and the measured D*~ 
momentum vector can be calculated from 4-momentum conservation in the B° — > D* + D*~ decay 

m% + m 2 D , - m 2 Dl - 2E B E D * 
cosiW = 9| M ; • (5) 

The B four-momentum is therefore determined up to the azimuthal angle around the D*~ direction. 
However, an arbitrary choice of this angle (e.g. cos 4>bd* = 0) introduces only a negligible spread 
(of the order of 1.5 MeV/c 2 ) in the missing mass distribution. 



4.2 Event Selection 

Random D*~ — 7 combinations are suppressed requiring | cos $bd* \ < 1-2. To reject events from 
continuum, we require the ratio of the second to the zeroth Fox- Wolfram moment (R2) .8 to be 
less than 0.3. 

Candidates for D*~ are reconstructed in the D°ir~ mode, using D° decays to K + ir~ , K + tt^-k + t:^ 
K + ir~ir°, and Kgir + Tr~ , here listed in order of decreasing purity. The x 2 probabilities of both the D° 
and D* vertex fits are required to be greater than 1%. The D*~ momentum in the T(AS) frame must 
satisfy 1.4GeV/c < pg^ s < 1.9GeV/c. Moreover, we require the reconstructed mass of the D° to 
be within 3 standard deviations <Jm D o of the nominal value m D \^ G jj, and Q B * = m B * —m D o —m^- 
to satisfy Q\ Q < Qd* < Qhi, where the limits Q\ Q = 4.10 to 5.20MeV/c 2 and Qhi = 6.80 to 
7.90MeV/c 2 are chosen around the nominal value Qf)P G = 5.851 MeV/c 2 depending on the D° 
decay mode. Kaon identification is required for modes K + tt~tt° and K + tt~tt + it~ . For mode 
Kg-K + TT~, the invariant mass of the tt + it~ from the K® decay is required to lie within 15MeV/c 2 
of the nominal mass and its flight length must be greater than 3 mm. 
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If more than one D*~ candidate is found, for each D decay mode we choose the best one based 
on the quality of the slow pion track and on the minimum value of x 2 = [(Qd* — Qd? G )/< j Q d *] 2 + 
[(m D o — nipo) I 'c m£)0 ] 2 , where cq d * is the measured resolution on Qd*- Finally, if candidates 
from different D° decay modes are present, we select the one with the best expected purity. 

Photon candidates are chosen from energy releases in the electromagnetic calorimeter not associ- 
ated with any charged tracks. In order to reduce random associations, we reject photon candidates 
which form, in combination with any other photon in the event, a 7r° whose invariant mass is 
between 115 and 155MeV/c 2 and whose momentum in the CM is greater than 200MeV/c. The 
selection of photon candidates is based on the optimization of the statistical significance of the 
observed signal (S/y/S + B, where S and B are the number of signal and background photons 
respectively), using Monte Carlo events. We require a minimum photon energy in the T(4S) CM 
fiCMS Q £ MeV, a minimum cluster lateral moment of 0.016, and a minimum Zernike moment 
|10j of order {2, 0} of 0.82. In about 10 % of the events, more than one photon is selected. In these 
occurrences we choose the one that maximizes the value of a likelihood ratio based on four photon 
variables (Ej, E^ MS , N cry , LAT), where Ey is the photon energy in the laboratory frame and N cry 
is the number of calorimeter crystals in the electromagnetic shower. 

4.3 Signal Yields 

The signal reconstruction efficiency is determined from a Monte Carlo sample of B° —* D* + D*~ 
events by performing an unbinned maximum likelihood fit to the missing mass distribution. The 
signal peak is well described by a Gaussian probability density function (p.d.f.), while the back- 
ground, which is mainly due to random D*~ — 7 combinations, is parametrized with the function 
-B(m-miss) = Bq(1 — e - ( mmiss ~ mmax )/ 6 )(m m i ss /m max ) c , where m max is the end point of the missing 
mass distribution. In the fit, we allow seven parameters to vary: Bq, b, c and m max in S(m m i ss ), 
and the mean, width and area of the signal Gaussian. We perform a single fit to all D° decay 
modes; the sum of the branching fraction-weighted efficiencies for the four reconstruction modes is 
computed from the number of signal events fitted in the range |m m i ss — rriD a \ < 45MeV/c 2 . The 
result is (eB) = (e< • Bf°) = (5.11 ± 0.03) x 10~ 3 . 

We have validated the fitting technique and the method of extracting the signal on the generic 
Monte Carlo sample. The distribution of the missing mass of partially reconstructed B° candidates 
is shown in Fig. for B°B° (including signal), B + B~, and continuum Monte Carlo events. From 
the signal yield, using Eq.EJwe obtain the result B(B° -> D*+D*-) = (1.962 ± 0.036) %, which is 
consistent with the value (1.97%) used in the generation of the Monte Carlo sample. 

Figure ^3 shows the missing mass distribution in our data sample. The same fitting procedure 
is used to extract the number of signal events. In the fit we allow all parameters to vary, except 
the width of the Gaussian signal, which is fixed to the value determined from fitting the signal 
Monte Carlo missing mass distribution. The result of the fit to the missing mass distribution is 
shown in Fig. \J]p. The signal yield is 7414 ± 345 events, corresponding to a branching fraction 
B{B° -> D*+D*~) = (1.854 ± 0.086) %. The result is stable over different run periods. 

5 FULL RECONSTRUCTION ANALYSIS 

5.1 Signal Extraction 

The full decay chain used for the B° — > D* + D*~ exclusive reconstruction is B° — > D* + D*~ — > 
(Df^)(D°7r~) with D° decaying in the four modes listed above, and Df — ► </>7r + — ► K~K + ir + . 
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Figure 1: Missing mass distributions in the Monte Carlo (a) and in the data sample (b). The 
B(m m i ss ) background fit (dashed line) and -B(m m i ss )+Gaussian total fit (solid line) are superim- 
posed. In (a) the different background components are also overlaid; starting from the bottom: 
uu-dd-ss, cc, B + B~ , B°B° (including signal). 

After applying selection cuts on the D* + and D*~ candidates, the combination with the smallest 
value of \AE\ = \{Eo* + Ed*) — -Ebeaml (ah quantities defined in the CM frame) is selected. Finally, 
the number of fully reconstructed B° is obtained from a fit to the spectrum of the energy-substituted 

mass mES = \A^beam ~~ (.Pd* + Pd*) 2 , where po* and po* are the D*~ and D* + momenta, again 
in the CM frame. 

5.2 Event Selection 

The selection of D*~ candidates, and most of the requirements on photon candidates in the full 
reconstruction analysis, are identical to those in the partial reconstruction. In the full reconstruction 
the background level is very small. We can therefore relax the requirement on the minimum 
photon energy in the center-of-mass system E^ MS , thus maximizing the statistical significance of 
our sample. 

(ft candidates are reconstructed from two charged tracks, with at least one track satisfying 
stringent kaon selection criteria. We use the helicity angle $k, defined as the angle between the 
kaon direction in the (ft rest frame and the (ft meson direction in the Df frame, to further suppress 
background. 

Df candidates are formed by combination with an additional track, with charge opposite to the 
slow pion from the D*~ decay. A mass window of ±50 MeV/c 2 around the nominal mass [Q is 
required. D*~ and D* + mass constraints are finally imposed in order to improve the AE resolution. 

At the end of an optimization procedure based on generic Monte Carlo which minimizes the 
overall statistical error and the peaking background contribution, we require the mD* s — ttid s mass 
difference to be between 0.125 and 0.160 GeV/c 2 , the reconstructed (ft mass to be between 1.0077 
and 1.0347 GeV/c 2 , |costf^| > 0.35, K? MS > 0.09 GeV, and \ AE\ < 0.05 GeV. 
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5.3 Signal Yields 



We determine the selection efficiency fitting the mES distribution of the signal Monte Carlo sample 
with a Crystal Ball function defining the number of signal events as the integral of this p.d.f. in 
the signal region 5.27 < hies < 5.29 GeV/c 2 . Summing the branching fraction-weighted efficiencies 
over the four D° reconstruction modes (see Sec. l4.2|) yields (e'B D °) = (6.28 ± 0.24) x 10~ 3 . 

Monte Carlo studies show the presence of a peaking contribution due to real B° — > D* + D*~ — > 
(Df'y)(D°7r~) events, where either the D° does not decay to the reconstructed modes, or the Df 
does not decay to 4>tt + . 

The total mES distribution is fitted with the sum of a Crystal Ball and a threshold function 
|12j . accounting for the combinatorial background. All parameters are allowed to vary in the fit, 
except the end point of the threshold function, which is fixed to 5.29GeV/c 2 . The signal yield is 
obtained from the integral of the Crystal Ball p.d.f. in the signal window, after subtraction of the 
peaking component discussed above. 

The fit procedure is first checked on generic Monte Carlo samples, where no bias is found in 
the full reconstruction analysis technique. 

Fig. [21 shows the fit to the total data sample and the corresponding yield. In order to subtract 
the peaking background from this yield, we scale the number of peaking events previously estimated 
from the Monte Carlo sample to the luminosity of the data sample. We then apply a correction 
factor to take into account the fact that the peaking events come from real B° — > D* s + D*~ decays 
and the measured value of this branching fraction in data (1.85%) is slightly lower than in the 
Monte Carlo. In addition, the peaking component coming from events with not decaying to 
<f)ir + , is rescaled using the measured (Df — > </>7r + ) branching ratio, with an iterative procedure. 
The resulting number of peaking background events expected in the data sample is 35, with a total 
uncertainty of 6 events which will be included in the systematic error. After subtraction of the 
peaking background events, the final yield on data is: 

Stot = (212 ± 19). (6) 

From equation H we therefore determine B 2 = B{B° -> D*+D*-) ■ B{D+ -> <pir + ) = (8.71 ± 
0.78 sta t) x 1CT 4 . The result is stable over different data-taking periods. 

6 SYSTEMATIC STUDIES 

The main sources of systematic uncertainties on the B° — > D* s + D*~ branching fraction measure- 
ment are listed in Tabled The Monte-Carlo-statistics uncertainty is due to the statistical error on 
the efficiency determination. The uncertainty due to the use of a fixed width for the signal Gaus- 
sian is estimated from the spread in the fit results allowing the width to vary, and simultaneously 
fixing in the fit the background shape as determined in several different ways. We conservatively 
assign an error of 4.5 %. The systematic uncertainty due to tracking efficiency is 0.9% per track 
and 1.6% for the soft pions from D*~ decays. The systematic error due to the isolated photon 
reconstruction efficiency and particle identification are evaluated using control samples. We find 
a 7% difference in the overall selection efficiency between our samples with complete longitudinal 
or transverse polarization in the B° — > D* + D*~ decay. The uncertainty due to the dependence 
on polarization is computed taking into account the experimental measurement of the fraction of 
longitudinal polarization, Tl/F = (51.9 ± 5.7)% |13j . Finally, the uncertainties on the D , D*~ 
and D* + branching fractions [T] are propagated throughout the analysis. 
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Figure 2: Total fit to the data sample. The Argus component parametrizing the background is also 
shown (dashed line). 



Table 1: Fractional systematic uncertainties 
surement. 



for the B° — ► D* + D* branching fraction mea- 



Source 


Error (%) 


Monte Carlo statistics 


0.6 


p.d.f. modelling 


4.5 


B counting 


1.1 


Tracking efficiency 


2.4 


Soft pion efficiency 


1.6 


Vertexing 


2.0 


Photon efficiency 


4.6 


Particle identification 


0.9 


Polarization uncertainty 


0.8 


D° branching fractions 


3.2 


B{D*- -> D°tt-) 


0.7 


B(D*+ - L»+ 7 ) 


2.7 


Total systematic error 


8.6 
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Table 2: Sources of systematic error (%) for the determination of B(DJ — > <pir + ). 





Source 


Error (%) 


Partial rec. 


Monte Carlo statistics 


0.6 




p.d.f. modelling 


4.5 


Full rec. 


Monte Carlo statistics 


3.2 




Tracking efficiency 


2.6 




Particle identification 


0.9 




Peaking background 


2.8 




Combinatorial background 


2.9 




B(<j> -► K~K+) 


1.2 




Total 


7.5 



Some systematic uncertainties, namely B counting, tracking efficiency, soft pion efficiency, pho- 
ton efficiency particle identification, polarization uncertainty, B(D*~ — > D°7r~) and B(D* S + — > 
Df^f)) cancel in the ratio (Eq. EJ). All other sources are listed in Table |2j The error on peaking 
background is due to the Monte Carlo statistics and to the uncertainty on the relevant D° and Df 
branching ratios; the error from the combinatorial background is estimated using the AE sideband 
data sample as an alternate way of computing the number of background events under the peak. 
The error on the B(4> — ► K~K + ) is taken from ^j. 

7 SUMMARY 

A measurement of the B° — ► D* + D*~ branching fraction is performed, using data corresponding 
to an integrated luminosity of 112.3 fb _1 , with a partial reconstruction technique. Including the 
systematic errors discussed in the previous section we obtain: 

B (B° D*+D*-) = (1.85 ± 0.09 stat ± 0.16 syst )%. (7) 

This preliminary result is compatible with, and improves on the precision of previously published 
experimental results ^EH> an d may be compared with the most recent theoretical results based 
on the factorization assumption 0: B(B° -> £>* + £>*~) t hcor = (2.4 ± 0.7)%. 
The preliminary Df — ► <j)7r + branching fraction result is 

B(D+ -> 0vr+) = (4.71 ± 0.47 stat ± 0.35 syst )%, (8) 

This new determination improves on published results [21 El 113 • 

8 ACKNOWLEDGMENTS 

We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the 
excellent luminosity and machine conditions that have made this work possible. The success of 
this project also relies critically on the expertise and dedication of the computing organizations 
that support BABAR. The collaborating institutions wish to thank SLAC for its support and the 
kind hospitality extended to them. This work is supported by the US Department of Energy and 



14 



National Science Foundation, the Natural Sciences and Engineering Research Council (Canada), 
Institute of High Energy Physics (China), the Commissariat a l'Energie Atomique and Institut 
National de Physique Nucleaire et de Physique des Particules (France) , the Bundesministerium fur 
Bildung und Forschung and Deutsche Forschungsgemeinschaft (Germany), the Istituto Nazionale 
di Fisica Nucleare (Italy) , the Foundation for Fundamental Research on Matter (The Netherlands) , 
the Research Council of Norway, the Ministry of Science and Technology of the Russian Federation, 
and the Particle Physics and Astronomy Research Council (United Kingdom). Individuals have 
received support from CONACyT (Mexico), the A. P. Sloan Foundation, the Research Corporation, 
and the Alexander von Humboldt Foundation. 

References 

[1] Particle Data Group, S. Eidelman et al, Phys. Lett. B592, 1 (2004). 

[2] The CLEO Collaboration, M. Artuso et al, Phys. Lett. B378, 364 (1996). 

[3] J. L. Rosner, Phys. Rev. D42, 3732 (1990). 

[4] The CLEO Collaboration, G. Brandenburg et al, Phys. Rev. Lett. 80, 2762 (1998). 
[5] The CLEO Collaboration, M. S. Alam et al, Phys. Rev. D50, 43 (1994). 
[6] Z. Luo and J. L. Rosner, Phys. Rev. D64, 094001 (2001). 

[7] The BABAR Collaboration, B. Aubert et al, Nucl. Instrum. Methods A479, 1 (2002). 

[8] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 

[9] A. Drescher et al, Nucl. Instrum. Methods A237, 464 (1985). 
[10] R. Sinkus and T. Voss, Nucl. Instrum. Methods A391, 360 (1997). 
[11] D. Antreasyan, Crystal Ball Note 321 (1983). 

[12] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B241, 278 (1990). 
[13] The BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. D67, 092003 (2003). 
[14] FNAL E687 Collaboration, P. L. Frabetti et al, Phys. Lett. B313 253 (1993). 
[15] BEPC BES Collaboration, J. Z. Bai et al, Phys. Rev. D52 3781 (1995). 



15 



